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Abstract 

We propose a novel, potentially useful generating technique for constructing 
exact solutions of inflationary scalar field cosmologies with non-trivial poten- 
tials. The generating scheme uses the so-called superpotential and is inspired 
by recent studies of similar equations in supergravity. Some exact solutions 
are derived, and the physical meaning of the superpotential in these models is 
clarified. 



Inflaton driven cosmological expansion, although not without some problems, still 
represents probably our best bet to describe the evolution of the early universe. On 
the other hand, the dynamics of the present day universe, we used to believe until 
now, is correctly described by a standard decelerating model. The new observational 
evidence [[I], |2[], however, points in favour of spatially flat weakly accelerating universe 
even at present. This weak acceleration must be induced by a matter with negative 
pressure. 

Some new models |3|, |J were proposed recently to account for the presence and for 
a possible present time dominance of energy density with negative pressure. Apart 
from usual radiation or matter, these models contain an additional ingredient either 
in the form of a constant vacuum stress, or quintessence, usually taken to be a scalar 
field slowly rolling down some non-trivial potential. In these model universes both 
the ordinary matter, the adiabatic perfect fluid, and the scalar field, are equally 
important, and one of the main enigmas to be resolved in this context is why the 
energy density of the additional quintessence matter almost co-incides with the usual 
matter density today |3], |5|. It follows, that the main player in these models becomes 
the shape of the scalar field potential, which due to scaling properties of the energy 
density of the scalar field || [7|, || determines much of the model phenomenology. 
In these scenarios one may divide the history of the universe into distinct periods 
depending on which source, the scalar field, or matter (radiation) dominates the 
expansion. Initially, to allow for the fast inflation, the universe should be dominated 
by the scalar field, then the radiation and matter energies take over, and then at 
present, the universe must start accelerating, signalling that the scalar field energy 
becomes important once again. 

During the periods of the scalar field domination, one may discard the other forms 
of matter and just concentrate on the behaviour of the scalar field driven expansion, 
therefore, the study of the scalar field models with different self-interaction terms 
becomes of great interest by its own right. Although many qualitative and numerical 
studies of scalar field models were undertaken during the last years due to the advent 
of inflationary cosmology, the number of exact solutions of Einstein field equations 
with self-interacting scalar fields, even in the presence of high symmetry, is rather 
poor. 

The main purpose of this Letter is to present a simple generating technique for 
solving coupled Einstein-scalar field equations with non-trivial potentials. We limit 
our discussion to the case of spatially flat isotropic universe, anticipating, that in 
more complex situations, either in presence of spatial curvature, or with extra matter 
fields, the method we employ doesn't seem to be too promising. 

To this end we consider spatially flat FRW line element given by: 





where A is a function of time alone. 

The Einstein Equations with the scalar field stress tensor are given by: 
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R»u = </>/A - 9^V{(f)), (2) 

and may be cast after some algebraic manipulations into the two following independent 
equations 

2A + <p 2 = (3) 
SA 2 - = V (4) 

The energy conservation equation is identically satisfied when (^) and (|]) hold. We 
have set the constants appearing in the Einstein equations to 1, and are using the 
usual GR normalisation for the scalar field as follows from (2). 

Given a potential V(4>) one may try to solve these coupled differential equations. 
This happens to be rather difficult a task. Some years ago Ellis and Madsen || have 
proposed a way of "solving" these equations by first assuming some particular form 
for the function A(t). Differentiating A(t) twice one may readily read from (|3|), 
and then you obtain V(t) from ([|). In other words, you just assume the expansion 
you wish, and then you read the potential after doing some calculus, "solving" the 
Einstein equations from "left to right". Since what one is really interested in is the 
potential as a function of you must invert <fi(t) into and then substitute this 
into V(t) to finally get V(<p(t)). This may be easily done for some simple choices of 
the scale factor a(t) = exp A(t), those like a(t) cx t n , or a(t) oc exp (kt), however if 
the scale factor is more involved the procedure of inversion of may become very 
complicated. The main flaw in this scheme, however, is that even if one manages to 
invert the above expressions, the potential V(<p) one finally finds may happen to be 
completely irrelevant. 

A different approach to solve this system of equations is suggested by their sim- 
ilarity to those ones arising in supergravity [|K|. Here, one introduces the so-called 
superpotential W(4>), such that 

V = -2(^) 2 + 3W 2 , (5) 

along with 

W{(j)) = A (6) 

and 

dW 

In 5-D gauged supergravity, the potentials of the form (||) occure naturally, but 
here, we will merely use W(<f)) as a solution generation function. It is straightforward 
to see that the Einstein equations are immediately satisfied by this choice, and more- 
over, the solution spaces of both sets of equations (|3|, [|) on one hand, and (|5], |B| and 
0), on the other, are identical. The strategy for solving the equations and |J) then 
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might be the following: Choose the superpotential W{4>). This immediately defines 
the shape of the potential, being the principle advantage of this generating scheme. 
Next, from (0) we can find t(<f>) and invert it to (j){t). The function A(t), on the other 
hand, is then obtained by the quadrature from @. The inversion is still a technical 
nuisance of the method, nevertheless, the transparency of the form of the potential 
is worth the effort. 

We now apply the above outlined scheme to obtain some exact solutions. First, we 
check that the superpotential W = cexp kef) which defines the exponential potential 

V(<f>) = c 2 (3 - 2A; 2 ) exp(2A;0) 



gives an expected power low scale factor JH 

i 

a = tvP, 

with the scalar field being 

(f) = \\\ogt + log2c£; 2 ] 
k 

Interestingly enough, a simple constant shift in W{4>) 

W = c exp kep + b, 

changes the potential to: 

V(4>) = c 2 (3 - 2A; 2 ) exp(2A;0) + 6c6exp (k<p) + 36 2 , 

with the scalar field being unchanged, but introducing a new behaviour into the scale 
factor 

a = exp (btftik 1 

It is as if we have added a cosmological constant to the potential and the scale factor, 
which had a power low behaviour before this addition, has been enhanced ( multiplied) 
by an exponential de Sitter term. 

We now consider a power low behaviour for W(4>): 

W(<t>) = k(j) n + b 

The scalar field potential is given then by: 

V{(j>) = 3k 2 <j) 2n - 2{nkfct) 2n - 2 + Qkbct) n + 36 2 , (8) 

and integrating ( and inverting) ([F|) we get for the scalar field ( for 

= [2kn{n - 2)t\*=* , (9) 
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Resolving (0) for the metric function we find 



1 n 2 

A = k(l- -n)[2kn{n - 2)}~t~ + bt 



In the case n = 2 we have 



and 



<p = exp (— 4kt) 
A = — exp(-8fct) +bt 



(10) 



(11) 



(12) 



Needless to say that the potential given by is quite rich and represents many 
interesting particular cases: for n — 1 and b = we have a typical case of quadratic 
potential, whereas for n = 2 and b = we have a quartic potential resulting in 
a double exponential scale factor a. Tracker type solutions may be obtained for 
negative powers of n < in W(<f>) @j. In terms of the superpotential, the scaling 
index £ of Ref. M is given by 



(13) 



here prime denotes derivative with respect to 0. It is easy to see that the scaling £ 
is constant only for a constant or for an exponential superpotential which results in 
exponential or power low behaviour of the scale factor respectively. Alternatively, if 
we write the scaling behaviour as oc l/a m , we then have that 



m 



WW 

w{<p) 



here one must assume that the "adiabatic" index of the perfect fluid model describing 
the scalar field is constant H. The identification of the perfect fluid and the scalar 
field, just to recall, is achieved by identifying the velocity potential of the fluid with 
the scalar field. 

The scalar field energy density can also be expressed in terms of the superpotential 
and simply becomes 



p=-tf + V( ( p) = 3W 2 ^), 



(14) 



hence, we arrive at quite an amusing result that the superpotential W(4>) has a rather 
simple physical meaning as the square root of the energy density of the scalar field. 

Phenomenologically, the cosmological models accelerate or decelerate depending 
on the sign of the second derivative of the scale factor. Expressing the acceleration 
conditions in terms of W gives accelerating models whenever 



W 2 



w{<p)_ 



>o, 
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otherwise we have a "standard" decelerating expansion, and it is not difficult to see 
what happens with each model. We have checked that depending on the power n and 
on the constants k and 6 one may construct with the potential (8) models to ones taste: 
ever inflating models, ever decelerating models, or models where the deceleration and 
inflation interchange several times. 

Our final example, just before closing, is the potential obtained from W = k cosh b(f>, 

V = 3k 2 cosh 2 (60) - 2{kb) 2 sinh 2 (60) (15) 

In this case the scalar field is 

= ^-Arctanh[exp (b 2 kt)), (16) 

and the metric function A may be integrated to give 

A= ^A;(coth 2 (6 2 f) + l)t, (17) 

and of course one may proceed further with generating variety of models. Therefore, it 
becomes relatively easy to construct exact solutions for different shapes of potentials. 
Moreover, since in principle one may express all the relevant dynamical information 
about scalar field models via the function W{4>) and its first derivative, this function 
may serve as an interesting dynamical variable to study the models analytically as 
well as doing some numerics, but this is away from the scope of this note. 

To sum up, we have presented simple generating procedure to construct exact 
solutions for scalar field isotropic and spatially flat cosmologies with potentials of 
different shapes. It is hoped that this method will become fruitful for generating and 
studying models of physical interest, as well as it is hoped that the superpotential 
function W((f>) might turn a useful tool for studying dynamical properties of these 
models. 

This research is supported by the University of the Basque Country Grant UPV 
122.310-EB150/98, General University Research Grant UPV172. 310-G02/99 and 
Spanish Science Ministry Grant PB96-0250. 
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